Theory for the couphng between longitudinal phonons and intrinsic Josephson 

oscillations in layered superconductors 
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In this publication a microscopic theory for the coupling 
of intrinsic Josephson oscillations in layered superconductors 
with longitudinal c-axis-phonons is developed. It is shown 
that the influence of lattice vibrations on the c-axis transport 
can be fully described by introducing an effective longitudinal 
dielectric function e^^iuj). Resonances in the /-IZ-character- 
istic appear at van Hove singularities of both acoustical and 
optical longitudinal phonon branches. This provides a natu- 
ral explanation of the recently discovered subgap structures in 
the /-V-characteristic of highly anisotropic cuprate supercon- 
ductors. The effect of the phonon dispersion on the damping 
of these resonances and the coupling of Josephson oscillations 
in different resistive junctions due to phonons are discussed 
in detail. 

74.80Dm, 74.50-Hr, 74.25Kc, 74.25Jb 



I. INTRODUCTION 

The c-axis transport in the highly anisotropic cuprate 
superconductors Tl2Ba2Ca2Cu30io+5 (TBCCO) and 
Bi2Sr2CaCu208+5 (BSCCO) can well be described by a 
stack of Josephson junctions between the superconduct- 
ing Cu02-multi-layers. This intrinsic Josephson effect 
can be seen in the multi-branch structure of the I-V char- 
acteristic, each branch corresponding to pa.. well-defined 
number of junctions in the resistive state.trB The intrin- 
sic Josephson effect is observed also in the behaviour of 
the material in extergial magnetic fields and under mi- 
crowave irradiationJ^I'B 

Recently subgap structures in the /-T^-characteris- 
tic have rbepn discovered as intrinsic properties of the 
matcrial,Lru which could be explained by the cou- 
pling between the intrinsic Josephson oscillations and 
phonons.ErQ This interaction is mediated by the oscillat- 
ing electric field in the Josephson junction, which excites 
vibrations of char|ged ions in the material. In our previ- 
ous investigationsBLI we used a simple model of a system 
of damped harmonic oscillators in order to describe the 
dynamics of oscillating ions in the barrier. We were able 
to derive an analytic expression for the dc-current den- 
sity j{V) as function of the dc- voltage V for one resistive 
junction 



3{V) = j,p(F) 



£2 



(1) 



where the voltage V is related to the Josephson oscilla- 
tion frequency lo hy V — TiLj/{2e) and e{Lo) = ei(w) -t- 
162 (ijj) is the dielectric function of the oscillating ions. 
From this result it can be seen that the /-F-curve has a 
maximum at the frequency (voltage) where the real part 
ei of the phonon-dielectric function vanishes, which cor- 
responds to a longitudinal eigenfrequency of the phonon 
system. 

With appropriate values for the Josephson plasma fre- 
qency Wp, the quasiparticle conductivity cr, the critical 
current denstiy jc, and the frequencies, dampings and 
oscillator strengths of phonons in the dielectric function 
we were able to fit the experimental results for the sub- 
gap structures in the /-y-curves perfectly. In addition 
to this, the identification of the maxima of the structure 
with phonon frequencies provides a natural explanation, 
why the position of the resonances is completely indepen- 
dent of temperature, magnetic field and the geometry of 
the probe. The voltage positions of these resonances £&«. 
also be related to structures in optical experimentd^'EHl 
in particular to the reflectivity for oblique incidencetJ . 

Despite the clear evidence of the effect a general the- 
ory for the coupling between Josephson oscillations and 
phonons is still missing and shall be presented in the fol- 
lowing. In particular, in Refs. |6||^ the phonon frequencies 
and the oscillator strengths in the dielectric function were 
not specified. The clarification of this point and the ex- 
tension to a general lattice dynamical model including 
acoustical and optical branches will be one of the main 
topics of the following investigation. 

Another topic concerns the coupling of Josephson os- 
cillations in different resistive junctions. We will show 
how phonons lead to a coupling of the phases of Joseph- 
son oscillations in different resistive junctions. This will 
be discussed in detail for two resistive junctions and gen- 
eral results will be given for large stacks of resistive junc- 
tions. A phase-locking in a stack of Josephson junctions 
is important for applications of such systems for high- 
frequency mixers and detectors. 

The excitation of phonons by Josephson oscillations 
in conventional single Josephsoii-junctions has been ob- 
served already a long time ago.t3 Also in the curves 
of break-junctions of cuprate-superconductorsEj struc- 
tures due to phonons might have been identified. The 
physical mechanism described here can also be applied 
to these cases but our formalism is particularly suited to 
treat stacks of Josephson junctions with phonons in the 
frequency range between the Josephson plasma frequency 



1 



and the gap-frequency. 

It is not our intention to explain the details of the 
/-^-characteristic of TBCCO and BSCCO by a reahs- 
tic lattice dynamical calculation. This is impossible at 
the moment for the complicated anisotropic supercon- 
ductors with variable doping showing this effect. Fur- 
thermore this would require a detailed theory of super- 
conductivity and the Josephson effect in two- and three- 
layer systems including information about the supercon- 
ducting bands, charge distribution, and charge suscep- 
tibility between the layers and inside the Cu02 planes 
which is not available yet. Therefore we discuss a simple 
model system with superconducting mono-layers where 
the conduction electron charge is distributed homoge- 
neously along the layers, and a lattice dynamical model 
with only two phonon bands showing already the basic 
features of a full theory which will be expected also for 
realistic systems. We start the discussion with a defi- 
nition of our model and a short derivation of the basic 
Josephson equations for a stack of Josephson junctions. 



II. JOSEPHSON EQUATIONS FOR A STACK OF 
JUNCTIONS 

We consider a system of TV superconducting layers sep- 
arated by insulating barrier material of thickness d form- 
ing a stack of Josephson junctions. We treat the super- 
conducting layers as homogeneous metal sheets with a 
uniform electron distribution along the layers. In this 
paper we treat only the case of a uniform tunneling cur- 
rent with a constant bias-current and neglect magnetic 
field effects due to the current flow. In this case all quan- 
tities are constant along the layers. Such an approxi- 
mation is reasonable for a stack of junctions which is 
short with respect to the Josephson penetration length 
but long enough to neglect finite size effects in the ionic 
polarisation. 

The tunneling current density j„ from layer n to n -I- 1 
creates (two-dimensional) charge density fluctuations 5pn 
on the layers related by the continuity equation 



jri(0 ~jn-l{t) = -Sp„{t). 



(2) 



These charge fluctuations create electric flelds E^{t) (in 
the c-direction) in the barrier between layer n and n -t- 1 
which are constant inside each barrier and are related to 
the charge fluctuations by the Maxwell equation: 



n' <n 



(3) 



(4) 



Assuming that the current density j„ for the first and 
last barrier is fixed by the bias current density j, then 



with help of Eqs. (H,§) the tunneling current density j„ 
in all the other junctions is related to the bias current 
density j by 



J ^Jn{t)+eomt)■ 



i5) 



The last term is the displacement current density related 
to the charge fluctuations on the layers. In the follow- 
ing we denote this term by Dn{t) ■= eoEP{t) in order to 
relate the present results to the usual notation of the re- 
sistively shunted Josephson junction (RSJ) model, how- 
ever, one should keep in mind the microscopic origin of 
this term. 

In the following we approximate the tunneling current 
by a superposition of a Josephson supercurrent density 
and a quasiparticle current density. Then we have for 
each junction: 



j = jcSin7„(t) +jqp{En{t)) + Dn{t). 



(6) 



The Josephson current density jc sinj„{t) depends on the 
gauge invariant phase difference 7„(i) between layers n 
and n -|- 1 at positions Zn and Zn+i- It is related to the 
average total electric field in the barrier 



En{t) 



Ez{z, t)dz 



by the second Josephson equation 



2ed 



j,,{t) = E„{t). 



(7) 



(8) 



Here small corrections to Eq. (||) for layered supercon- 
ductors which are discussed in Refs. p^ljl^ are neglected. 
For the quasi-particle current density we will use in the 
following an ohmic form jqp — aEn with a constant con^ 
ductivity a. The generalisation to more realistic formsQ 
is straightforward. 

The crucial point where the phonons come into play 
is the relation between the field Dn = ^qE^, which is 
created by the charge fluctuations on the superconduct- 
ing layers alone, and the average electric field Eq. (|^) 
En = E^ + En'^ which contains also the averaged field 
created by the ionic displacements in the barrier. 
This will be discussed in detail below. 

Before we do this let us summarize the most impor- 
tant parameters which characterize the Josephson sys- 
tem: The first one is the (bare) Josephson plasma fre- 
quency ujp defined by 



2edjc 
heo 



(9) 



The second one is the so-called characteristic frequency 
defined by 



2eVr 



(10) 
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Here Vc is the voltage where the quasiparticle current 
density equals the value jc- It is of the order of the su- 
perconducting energy gap and is a measure of the dissi- 
pative properties of the junction. In our simple model 
with a constant conductivity we have tOc — 2edjc/(ha). 
The ratio /3c — ^c/'^p is the McCumber parameter, 
which is /3c 3> 1 for the strongly anisotropic cuprate- 
superconductors. Moreover, for these materials there ex- 
ist phonons with frequencies in the range ujp < Wphon ^ 

Typically the time-dependence of the phase difference 
7(t) can be written in the so-called resistive state as 

-/{t) = 9 + Lot + S-/{t) (11) 

where uj = (7) is the time-average of the phase velocity 
which is non-zero for a junction in the resistive state. It 
determines the dc- voltage V = {E)d = h{'j)/{2e) across 
the junction. In the asymptotically stable state and for 
large values of the McCumber parameter /3c the oscillat- 
ing part Sj{t) is small and oscillates with the frequency 

LU. 



III. EXCITATION OF PHONONS BY 
JOSEPHSON OSCILLATIONS 



E D^f,{£,)ep{K'\qX) = uj^iqX)e^{K\qX) (14) 
defined by 

The force-constant matrix contains the quantum- 
mechanical short-range interactions but also the short- 
range and long-range Coulomb interactions (the latter 
being of the form ^ qaqp/q^) between the ions, but not 
the fields set-up by the conduction electrons on the layers. 
The eigenfrequencies u){qX) are therefore by construction 
the hare phonon frequencies in the absence of the con- 
duction electron charge fluctuations 6pn(t). They include 
possible renormalisation by intra-atomic electronic polar- 
isation. 

For the further discussion it is convenient to label 
the lattice dynamical unit cells by I — {lx,ly,lz) with 
Iz = n denoting the superconducting layer in which 
the lattice cell is contained (see Fig. |l|). Then the z- 
component of the position vector i?(^) for I — {lx,ly,n) 
does not depend on Z|| = {lx,ly) and we may write 

RzC^) = Rzi") = Rz{n) + R^{k). Furthermore the ori- 
gin of the unit-cell may be chosen on the superconducting 
layer, so that Rz{n) ~ Zn- 



Now let us turn to the discussion of lattice vibrations. 
Quite generally the lattice displacement of an ion of type 
K with mass M^, charge Z^, in unit cell I is determined 
by the following equation of motion 



Z^EP^Cjt)- (12) 



I'k' 



Here £'^(^|i) is the local dr iving field at the equilibrium 

position R{[) = R(l) + R{k) of the ion generated by the 
charge fluctuations Spn(t) on the superconducting layers. 
Note that the vibrating ions may be both in the barrier 
material and on the superconducting layers. The super- 
conducting electrons are assumed to move together with 
the ions of the layers. 

While the RSJ-equations for the phases are highly non- 
linear, the relations between lattice displacements and 
electric fields are linear, consequently we may analyse 
the response for each frequency ui separately. With a 
harmonic ansatz of the form Ua{'^\t) ~ Ua{l_)e~ 
obtain for the amplitude 



we 



ea{K\qX)e*p{K'\qX) 



(13) 



N ^ UJ^(qX)-LU^ 

Here w^(gA) and e{K\qX) are the eigenvalues and eigen- 
vectors of the dynamical matrix 



n+1 



l=(l„,n) 



FIG. 1. Labeling of unit cells. 



n 



n-1 



In our model for the superconductivity in the layers we 
have to assume that the mobile electronic charge is spread 
out uniformly along the layers, because a microscopic 
theory connecting the superconducting bands and the 
atomic sites in the Cu02-layers is still missing. There- 
fore the field E^{1_) does not depend on the position 



along the layer and has only a z-component. In Eq. (|l£ 
only modes with gn — appear and we may write for the 
displacement amplitude in z-direction of an ion of type k 
in barrier n: 

ez{K\qzX)e*{K'\qzX) 
27 — 5 (16 



Now let us specify the driving field i?f (") for the ions 
in more detail: Summing up the electric fields created by 
the charge fluctuations on the different layers we find: 
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J2 ^Pn' - J2 ^Pn') for K S barrier 



n'<n n'>n 



2i^( Z ^Pn' - J2 ^Pn') for >^ € layer. 



(17) 



Note that for ions inside the barrier the driving field does 
not depend on the position Rz (k) of the ion and is equal 
to the constant field introduced in Eq. (||): 

j7P/n \ _ i En for K e barrier 

I J - I 1 (^p + EP^_^) for K e layer. ^^^> 

In order to find a relation between the lattice dis- 
placements and the driving field we introduce the Fourier 
transformation 



Then 



with the g^-dependent effective charge 

~ , , r for k G barrier 



and consequently 



(19) 



(20) 



(21) 



(22) 



E'{q.). 



Thereby the vanishing of Z^ for qzd = tt and for ions 
on the layers refiects the fact that for alternating elec- 
tric fields in the barriers no net force is acting on the 
superconducting layers. 

Next we want to relate the driving field to the av- 
erage total electric field i?„ in the barrier because this 
field is connected with the phase difference 7„ by the sec- 
ond Josephson-equation (^). The microscopic total elec- 
tric field Ez{x,z) is the sum of the fields set-up by both 
the charge-fluctuations and the displaced ionic charges: 
Ez{x, z) = EP{x, z) + El°''{x, z). As EP{x, z) is assumed 
to be constant inside the barrier and independent of x 
we may replace the other two fields also by their aver- 
ages across the barrier and over one unit cell along the 
layer. For the averaged fields we have then the relation: 



E^, — En 



Ell 



(23) 



The field in the barrier can be calculated from the 
difference of the scalar potentials on layers n and n + 1 
produced by the ionic displacements averaged over the 
area of one unit-cell as 



Ell 



— ( E z,m:) 

CO^^c ^ KSbarrior 



(24) 



+1 E z^uzO+uzrt')))- 



At G layer 

Here Vc is the volume of one unit cell. The factor 1/2 
results from the fact that displacements of ions (with 
Qll = 0) on a layer do not contribute to the potential on 
the same layer. Note that this field for q|| — does not 
depend on the ionic displacements in other barriers and 
is closely related to the ionic polarisation in the same 
barrier. Defining a generalised polarisation by 

Elr =■■ -Pn/eo, (25) 

we may write for the displacement in the barrier 

Dn-.^eoEP^eoEn + Pn, (26) 

which has the usual form as in the macroscopic Maxwell 
theory. 

Going over to a Fourier transformation the relation 



Piqz)=xiqz,io)eoEPiqz) 



(27) 



between the polarisation and the driving field is obtained 
with 



xiQz,'^) = E 



Lu'^iqzX) - up- 



(28) 



and the oscillator strength 

"'^ e2(K|q^A)e:(K'|g^A) ~ 



\^{qz\)f = — Y^Z^{qz 



(29) 



The special combination of phase factors contained in 
■Z^kX^z) Eq. ( |2l| ) are a consequence of the different contri- 
bution of ions on and between the superconducting layers 
to the electric field in the barrier. 

Using e^EP{qz) — e^E{qz) + P{qz) we can solve for 
^fe): 



P{q. 



x{Qz,uj) 

1 - x{qz,uj) 



£oE{qz 



(30) 



The relation D{qz) — eQ€^^^{qz,uj)E{qz) defines an effec- 
tive longitudinal dielectric function 



1 - x{qz,uj) 



(31) 



This function has zeros at the eigenfrequencies w(qzA) of 
the dynamical matrix. Due to the form of the oscillator 
strengths Eq. (|2^) only modes with polarization in the c- 
direction contribute. As the electric field En is constant 
along the layers we have gj| = and only longitudinal 
modes in the c-direction couple, therefore the zeros of 
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£pi^{qz, ^) are exactly at the longitudinal eigenfrequencies 
of the dynamical matrix. 

In the case of a single dispersionless phonon-mode with 
frequency cul the function £py^{i^j) can be directly com- 
pared with the dielectric function used in Ref. ^. In fact, 
in this case ep^^{u!) can be written as 



(32) 



The form of the longitudinal dielectric function 



epf^{qz,uj)) Eq. (|3l|), which we have in- 



troduced here, is different from the transverse dielec- 
tric function e'^^{qx,u}) = e^^(fe,t^)) = 1 + x(9x,w) In 
both functions different eigenfrequencies and oscillator 
strengths enter, however in the limit q^ — ^ 0,qx ~^ 0, 
which is only relevant in optical experiments, the values 
of both functions are equal. 

Finally a comparison of our theory with theoretical in- 
vestigations in Ref. |l^ are in order. In principle there 
are two different electron-phonon coupling mechanisms, 
which may couple Josephson oscillations and phonons: 1. 
the electromagnetic interaction between the ionic charges 
and the charges of conduction electrons, 2. the depen- 
dence of the tunneling matrix element on lattice dis- 
placements. The first mechanism is considered in our 
work, the second in Ref. ^ Both mechanisms require 
a different theoretical treatment (on a diagrammatical 
basis the two mechanisms would correspond to differ- 
ent diagrams). It has been argued in Ref. |l^ that in 
the layered cuprate superconductors the charges of the 
ions in the insulating barrier between superconducting 
layers are unscreened and therefore have a strong in- 
teraction with conduction electrons in the Cu02-layers. 
We therefore considered this mechanism for our treat- 
ment of the coupling between Josephson oscillations and 
phonons. Though we did not write down a Hamiltonian 
for the interacting system our method nevertheless is a 
full microscopic theory which treats the electron-phonon 
interaction on a random-phase-type level by describing 
the interaction with the help of internal fields. This ap- 
proximation is sufficient as long as we do not want to 
consider the electron-phonon interaction inside the su- 
perconducting layers and treat exchange effects between 
different superconducting layers. 



IV. INFLUENCE OF PHONONS ON THE 
IV-CHARACTERISTIC 

According to the RSJ-like model derived in Eq. (^ the 
current density in junction n is 

j = jc sin7„(i) + aE„{t) + b^{t) (33) 

where Dn{t) = eoEP{t) is the electric field in junc- 
tion n set-up by the charge fluctuations of conduction 
electrons. As pointed out in the previous section this 



field can be expressed by the average electric field in 
the barrier and the generalized polarisation Eq. (25) as 
^ri(^) = ^oEn{t) + Pn{t). The polarisation has to be 
calculated self-consistently from the ionic displacements 
and depends linearly on the electric field. 

Let us discuss first the case of one resistive junction 
at n = in the middle of a large stack while all other 
junctions n ^ are in the superconducting state. Then 
as mentioned previously all the oscillations are governed 
by one frequency lo, and we can write for the phase for 
n = 0: 



while for n ^ we have 

7n(i) = On+S-f„{t). 



(34) 



(35) 



In the stationary state 6'yn{t) oscillates with the same 
frequency lu, 



S^nit) ^ dj„e~"^' + c.c. 



(36) 



Higher harmonics can be neglected for /3c ^ 1. In this 
limit the fluctuations 5^n{t) are small and we may use 
the expansion 

sin 70 (i) ~ sin(6'o + oJt) + cos(6'o + u}t)5-io{t), (37) 

while for n ^ Q we have 

sin7„(t) ~ sin6'„ -I- cos 6'„57„(t). (38) 

The bias current density j on the l.h.s. of the RSJ- 
equations Eq. ( ^3[ ) is time-independent and equal for all 
junctions, while the quantities on the r.h.s have both 
time-independent and oscillating components. 

Let us discuss the equations for the non-resistive junc- 
tions [n 7^ 0) first. Here the dc-component is: 



J = ]c sm ( 



(39) 



This fixes the constant part of the phases in the non- 
resistive junctions and relates it to the bias current. 
For the oscillating part of Eq. (B3) one obtains: 



= jcCos6i„57„(i) -f a— <57„(i) -f I)„(i) (40) 
Zed 



or 



2ed ■ 



= Lu'Sjnit) + -djn{t) + —D„{t) (41) 
with the reduced Josephson plasma frequency 




(42) 



Now we discuss the resistive junction at n = 0. Keep- 
ing only the lowest harmonic we find 
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sin 70 W - sin(0o + (^t) + ^{S-foe'^°). 



(43) 



The dc-component of the RSJ-equation Eq. ( p3D is there- 
fore given by 



with 



where V is the dc-voltage of the resistive junction and 
jqp{V) = uEdc is the quasiparticle current density. 
For the oscillating part one finds: 



= jc sm{9o + Lot) + a—S^oit) + Do{t). (45) 
Zed 

The two equations (^, ^5|) can be combined to one 
inhomogeneous linear differential equation for all n 

LulSjnit) + -S^nit) + ^Dnit) = f„{t) (46) 
£0 "Co 

with 

f (i\ - f ^l^loit) - i^l sin(e'o + Ljt) for n = :.rj. 
J^'^^>-\0 forn^O. 

Assuming a time dependence of the form e^*'^* for all 
oscillating quantities we have 

ColSjn + ^^Sjn + ^{-iuj)Dn^ fn (48) 

with 



= h^p'^To-^e-^^" forn = (49) 
i for n 7^ 0. 

In order to incorporate the non-local dependence of 
the polarisation on the electric fields in different barriers 
a spatial Fourier representation of the form 



is introduced. Using the relation 



(50) 



D{qz) = eoepi,{qz,uj)E{q,) = ^Cphl 



(51) 

in Eq. yields 

G-\qz,uh{qz) ^ fo (52) 

with 



G ^{qz,uj) =ujl-iu— -uj'^epi^{qz,uj). (53) 



For the phase oscillation in the resistive junction follows 
^70 = X! " 9{^)fo (54) 



(44) Solving for ^70 



1 ^i^Jt, a 

2g \uj)-uj^ 



we obtain 



2 g ^{uj)-ujI 



(55) 



(56) 



(57) 



From this finally the following expression for the dc- 
current density as function of the dc voltage is obtained: 



2 w2 g(t^) 

3c h (^) 



(58) 



2 w2 + ii{Loy 

Here e(w) is a modified dielectric function 



where 



eph(t^) = J ^{lo) + ■ 



and 



(59) 



(60) 



(61) 




3.0 4.0 

(0/(271) [THz] 

FIG. 2. Analytical I-V-cuive for one resistive junction 
with one phonon resonance a.t lu — ujl- The arrows mark 
the hysteretic jumps found in current-biased experiments and 
numerical simulations. 
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This expression describes the de-current density as 
function of the dc- voltage V = Uuj/{2e). It has a maxi- 
mum for frequencies uj where the real part e(a;) vanishes, 
i.e. for longitudinal phonon frequencies. 

This can easily be seen if we consider the special case 
of one phonon mode without dispersion. Then eph(ti') = 



e(z,fc,w) 



9. ^n\v\1z^^) - 



9^ '^ph 



(65) 



i.k 



(cj), which is of the form Eq. (|3|). The corresponding The terms in brackets in Eqs. (|64|, |65|) are understood 



. . . M . . 

I-V-CMTYe is shown in Fig. which is calculated with the 

dielectric function 



1 



UJ, 



(62) 



Here an additional damping p of the phonon has been 
introduced in Eq. (^2|) in order to simulate the energy 
transfer into other junctions due the coupling of the ions. 
A peak in the /-T^-curve appears a,i uj = lo^. The width of 
the peak is determined by this damping together with the 
quasiparticle conductivity. The deviation from the quasi- 
particle current density vanishes at the pole of epj^(a;) at 

Wo = — f^-^. The rise at low voltages is due to the 

plasma resonance. For /?c ^ 1 the minimum of the I- 

y-curve is at ojret — (3/2)"'"^'*ijJp/y^e^^^^tijp). In current- 
biased experiments and corresponding numerical simula- 
tions the parts with negative differential conductivity are 
skipped hysteretically as indicated in the figure. 

Note that the denominator in the function J(w) Eq. 
( |6l| ) is the total q^- and cj-dependent longitudinal di- 
electric function of the coupled system of phonons and 
conduction electrons, 



{qz,^) 



iOZ, 



t^ijj' 



(63) 



Zeros of the real part of this function describe longitudi- 
nal collective modes in the system. On the other hand 
the resonances in the l-V curve appear at the hare longi- 
tudinal phonon frequency in the case of a narrow phonon 
band. The summation over in Eq. ( |6l| ) leads to an ef- 
fective damping of the resonances which is proportional 
to the phonon dispersion. The physical origin is the loss 
of energy by phonons from the resistive junction to the 
neighboring junctions. 

The result for the current-voltage characteristic can be 
generalised to the case of several junctions being in the 
resistive state, if we assume that all junctions oscillate 
with the same frequency lo. Denoting the subset of in- 
dices for the resistive junctions by / then for i (z I we 
obtain (for a derivation see the appendix): 



(64) 



P 



fee/ 



ell, 



eouj 



Si,k 



The dielectric function e(i, k,uj) is defined by 



as matrix inversions. The dc-voltage V is obtained, if one 
multiplies huj/{2e) by the number of resistive junctions. 
Note that the r.h.s. of Eq. ( |64| ) depends on the layer 
index i G /, while the l.h.s. is equal for each junction. 
From this equality the (relative) phases 9i in the different 
junctions can in principle be determined, which in turn 
provides an analytical expression for the I-V-cmve. 

In the case of two resistive junctions exactly two so- 
lutions exist with 6i — 9j and 9i = 9j + tt, respectively. 
In the general case several different solutions are found. 
The stability of these solutions will be checked by a com- 
parison with a direct numerical integration of the cou- 
pled equations of motions in the following section. It 
turns out, that those analytical solutions are most stable 
where the phases Oi of the oscillating Josephson junctions 
fit best to the pattern of lattice vibrations with the given 
frequency to. 



V. A SIMPLE EXAMPLE 

The theory in the preceding sections is developed for 
general lattice dynamical models and is valid within the 
assumptions we have made for the superconducting prop- 
erties: we treat only single-layer systems assuming a ho- 
mogeneous conduction electron charge distribution along 
the layers. An extension of the theory to more realistic 
systems is in principle possible, but requires to introduce 
the charge susceptibility of conduction electrons in the 
superconducting state and a generalisation of the Joseph- 
son theory to multilayer systems. But for this more de- 
tails of the electronic properties are required than are 
currently known about these materials. In addition to 
this, it is not possible to compare theoretical results for 
the /-^-curves with experiments in detail, as reliable lat- 
tice dynamical calculations for the strongly anisotropic 
systems BSCCO and TBCCO with two and three layers 
and variable oxigen content are not yet available. 

In the following we consequently only want to illustrate 
the main features of our theory in a simple toy model, 
which reflects some basic aspects of the real system. One 
of the main lattice dynamical property of these systems 
is certainly the existence of a longitudinal acoustical and 
(several) flat longitudinal optical bands which result from 
movements of groups of ions in the bacrijei; against ions 
in the Cu02-planes in the c-direction.llaE3 Such modes 
we simulate by the most simple lattice dynamical model 
consisting of two kinds of ions with ionic charges Zi, 
Z\j = —Zi and masses Mi, Mj,. The first kind {k — I) 
is placed on the superconducting layers, the second kind 
(k = b) in the middle of the barrier. The motion of ions 
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in the c-direction which is assumed to be uniform along 
the layers is approximated by a two-atomic chain-model 
with nearest-neighbor interactions in the c-direction: 

M,u{l) - fiuCY) + <D - Ml)) = ZtE'^il). (66) 

By choosing the masses very different a narrow optical 
band can be simulated. From a diagonalisation of the dy- 
namical matrix given by Eq. ( |66| ) the well-known eigen- 
frequencies uj^qzX) of the two-atomic chain are obtained. 
With help of the eigenvectors the oscillator strengths de- 
fined in Eq. ( p9|) are calculated, which are needed for the 
longitudinal dielectric function Eq. (^). 

The driving field on the r.h.s. of Eq. (|66| ) is the field 
set-up by the conduction electron charges on the super- 
conducting layers, which can be expressed by the (con- 
stant) field EP in the barrier between layers n and n + 1 
with the help of Eq. (|T^). The latter can be expressed 
by the average total electric field En in the barrier: 



where the polarisation Eq. (Ea) is given by: 



Pn 



^fzMi) + U{ucn+ucT)) 



(67) 



(68) 



These equations for the motion of lattice displacements 
have to be supplemented by the extended RSJ-equations 
Eq. (1): 

j = jcsin7„(i) + a-^7„(<) + to^ln{t) + Pn{t). (69) 
zed 2ed 

This model is used to calculate the /-^-characteristics 
in two ways: 1. The J-y-curves are calculated analyti- 
cally using the results Eq. (^) obtained by the Green's 
function method. Thus the peaks due to the phonon res- 
onances are obtained. 2. The coupled set of RSJ equa- 
tions Eq. (69) and phonon equations Eqs. (^, |8|) are 
integrated numerically by a Runge-Kutta method for a 
finite stack of Josephson junctions. Changing the bias- 
current gradually allows to follow the /-l^-curves as in 
the current-biased experimental situation and to repro- 
duce the hysteretic behaviour in particular. 

We start with the discussion of the first branch of the 
/-y-curve, where one junction is in the resistive state. 



the barrier contribute to the oscillator strength, the ions 
on the superconducting layers are inactive due to the 
factor 1 -I- exp(i(7^o?) in Eq. (pTf). These features will be 
illustrated in the following. 
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FIG. 3. Dispersion and oscillator strengths for the 
two-atomic chain model. Shown are the oscillator strengths 
A)| for the acoustical branch (dashed curve) and the op- 
tical branch (solid curve) and for the two cases of the heavy 
ions on the superconducting layer (M; > M^) and in the bar- 
rier [Ml < Mb). 



A. One resistive junction 

Quite generally the /-V"-curve is expected to have 
peaks at the van Hove singularities of the phonon disper- 
sion. Details, however, depend on the oscillator strength 
defined by Eq. ( p9| ) which enters the longitudinal dielec- 
tric function Eq. (pT]). In particular at the edge of the 
Brillouin zone for = 7r/d only the motion of ions within 



For the lattice dynamical model introduced above at 
Iz = /d only one type of particles is moving due to sym- 
metry: In the acoustical branch the heavier ion, in the 
optical branch the lighter ion is moving. If the heavier 
ion is on the superconducting layers (M; > Mf,) the oscil- 
lator strength vanishes at the end of the acoustic branch 
(see Fig. ^), and peaks are expected to appear in the 
/-F-curve at the two van-Hove singularities of the opti- 
cal branch. On the other hand, if the lighter ion is on 
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the superconducting layers (Mi < Mj,) then the oscilla- 
tor strength vanishes at the end of the optical branch, 
and peaks are expected at = tt/cJ from the acoustical 
branch and at = from the optical branch. 
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FIG. 4. /-^/-curves for one resistive junction with sub- 
gap-structures due to acoustical and optical phonons. 



phonon dispersion is fixed by the values of u}^Q{qz = 
0) = 2f{l/Mi + l/A'h) and the mass ratio Mi /Mb. 
A measure for the oscillator strength is the quantity 
fi^j, :— Z^/ {Mi^b^QVc)- In Fig. ^ we have used the follow- 
ing parameters: (ic = w^/wj = 375, Mi /Mb = 0.2 and 
Mb /Ml =0.2 respectively. The phonon frequencies are 
normalised to ojc and are given by ix'loI'Zz = 0) = 0.28wc. 
The oscillator strength is given by uJc/^f = 200 for 
Ml < Mb and uj^/^l = 200 for Mi > Mb. The corre- 
sponding phonon dispersion and the oscillator strengths 
are shown in Fig. |[ The values of the parameters (3c, 
i^Lo/^c are adapted to TBCCO, also the chosen value 
|f2p/wo = 0.08 for the oscillator streiigth compares well 
to the experimental estimate « O.lSllJ. The mass-ratio 
is chosen to obtain a sufficiently flat optical branch (cf. 
the band structure in Ref. ^^ ). 

For Ml < Mb the analytical results for the I-V-cmve in 
Fig. ^ show phonon peaks at the van Hove singularity of 
the acoustic branch at Qz — n/d and at the van Hove sin- 
gularity at = of the optical branch, the resonance at 
the van Hove singularity at Qz = n/d of the optical branch 
is suppressed. For Mi > Mb only structures due to the 
two van Hove singularities of the optical branch appear. 
In both cases the increase of the J-F-curve at low fre- 
quencies indicates the Josephsoriplasma frequency. The 
numerical results shown in Fig. ^ for Mi < Mb (here Wc 
= 13.1 THz) follow the analytical results in the regions 
of positive differential resistance perfectly, and otherwise 
show the hysteretic behaviour as seen in experiments. 
At low values of j/jc the I-V-cmve switches back to the 
superconducing state of the junction. 




0.0 I — ■ — ■ — ■ — ■ — ' — ■ — ■ — ■ — ■ — ' — ■ — ■ — ■ — ■ — ' — ■ — ■ — ■ — ■ — ' — ■ — -J 
0.0 1.0 2.0 3.0 4.0 

V [THz] 

FIG. 5. Comparison between analytical (solid line) and 
numerical results (rhombs) for the I -V -curve of one resistive 
junction with subgap-structures due to acoustical and optical 
phonons for A/; < Mb. 

This is illustrated in Fig. ^ where we have plotted re- 
sults for the I-V-cmve of the first branch in the two 
cases Ml < Mb and Mi > Mb- In our model the 



B. Two resistive junctions 

Another important effect of the coupling between 
Josephson oscillations and phonons is the synchronisa- 
tion (phase-locking) of Josephson oscillations in differ- 
ent resistive junctions, which would be absent without 
phonons in short junctions, which are homogeneous par- 
allel to the layers, we want to illustrate this for the case 
of two resistive junctions coupled by one narrow optical 
phonon branch. 

The second branch of the I-V-cmve which corresponds 
to two resisitve junctions has a rather complex structure 
already for one phonon band. This is shown schemati- 
cally in Fig. H If we denote the two dynamical states of 
the first branch by a and 6, then on the second branch 
with two resistive junctions either both junctions are in 
state a (label aa), both junctions are in state b (label bb) 
or one junction is in state a while the other junction is 
in state b (label ab). Note that in the latter case the os- 
cillation frequencies of the two junctions are different. In 
the case of well separated resistive junctions the voltages 
of a structure for a given bias current are determined by 
LUaa = '2uJa, ^bb = 2cjf,, uJab = cJq + cjf,. This is no longer 
true, if the resistive junctions are close to each other and 
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interact via phonons. Then the voltages in the second zone. The numerical results follow the in-phase solution, 
branch are slightly lower. Fig. |^c shows results for j — i + 3. 



1 . branch 



2. branch 




V 

FIG. 6. Schematic plot of the first and second branch of 
the J-V-curve with subgap structures due to one phonon. 

In the case of two resistive junctions i and j two solu- 
tions of Eq. (^) exist with phase differences 6i ~ 9j and 
Oi = 9j + TT, corresponding to in-phase and out-of-phase 
Josephson oscillations, respectively. Inserting these re- 
sults in Eq. two different J-T^-curves can be calcu- 
lated. Note that this formula applies only for the states 
aa and bb, because in the derivation we have assumed 
that the two junctions oscillate with the same frequency. 
A similar formula can also be derived for the state ab. 
In that case the phases are meaningless because the two 
junctions oscillate with different frequencies and are es- 
sentially decoupled. 

In Fig. ^ we show examples calculated for a narrow 
optical band with Mi /Mb = 10 and the light ion oscillat- 
ing in the barrier. Here the different analytical solutions 
are shown together with numerical results for neighbor- 
ing resistive junctions, j = i + 1 {a), and two resistive 
junctions separated by 1 or 2 superconducting junctions, 
j = z + 2 (b), j=z + 3 (c). 

It is plausible that in the case of neighboring resistive 
junctions out-of-phase Josephson oscillations {9i — 9j + 
tt) favor a coupling to phonons at the edge Qz = tt/c? of the 
Brillouin zone, while the coupling of in-phase oscillations 
{9i = 9j) is strongest for zone-center phonons at Qz = 0. 
This is shown in Fig. where the I-V-curve for the in- 
phase solution shows a peak at uj{qz — 0), while for the 
out-of-phase solution the current-maximum is at Lu{qz — 

TT/d). 

It can be seen that the numerical results in the dynam- 
ical state aa follow one of the analytical solutions before a 
hysteretic switch into state bb occurs (outside the figure) . 
This is verified in Fig. |^a where the numerical I-V-curve 
follows the analytical curve for 9i = S^+i -I- tt. In Fig. 
^ the in-phase solution with 9i = 9i^2 has maxima for 
voltages corresponding to frequencies of optical phonons 
at both Qz = n/d and qz = 0, while the out-of-phase solu- 
tion has a broad maximum in the middle of the Brillouin 
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FIG. 7. I-V-cuTve for two resistive junctions separated 
by (a), 1 (b), 2 (c) superconducting junctions. Shown are 
analytical /-^/-curves calculated for a narrow optical band 
for the in-phase, 9j — 9i — (solid line) and the out-of-phase, 
Oj — di = -K (dashed line) solution together with numerical 
results (dots). 



C. Several resistive junctions 

The results obtained for two resistive junctions show 
that the peak position in the I-V-cuive depends only 
slightly on the distance between the resistive junctions. 
More important is the fact that phonons are able to syn- 
chronise the phases of Josephson oscillations in different 
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resistive junctions. This is important for the use of such 
systems in high-frequency apphcations. We have checked 
this numerically for the case of a block of many resistive 
junctions. For frequencies lo close to a phonon cigenfre- 
quency uj{qz = Tr/d) at the edge of the BriouUin zone the 
Josephson field oscillations in neighboring junctions differ 
by a value close to tt. For frequencies close to u){qz — 0) 
near the F-point the Josephson oscillations are nearly in 
phase. These phase-locked dynamical states are reached 
from arbitrary initial conditions for the phases. There is 
also a synchronisation of Josephson oscillations for fre- 
quencies far away from phonon resonances. We did not 
yet investigate the stability of these states systematically 
but we expect these to be less stable than at frequencies 
close to a phonon resonance. 

VI. EXPERIMENTAL RESULTS 

Recently the explanation of the subgap resonances in 
Refs. with the phonon coupling mechanism pre- 

sented here could be well ccaifamcd by Raman measure- 
ments on the same samplesoEil and .infrared reflectivity 
experiments with grazing incidence, EilO where the lat- 
ter allows to determine both longitudinal and transverse 
modes (see Table ^. Small deviations of the order of 
~ 5 — 10% may be attributed to the fact that in optical 
experiments and in the intrinsic Josephson effect differ- 
ent averages over q of the dielectric functions are relevant. 
Note that in our theory modes which are Raman active 
at = may couple to intrinsic JosephsoB-oscillations 
also for qz ^ 0. Earlier experimental data,Elll3 which are 
obtained from polycrystalline samples, show the same 
qualitative behaviour, but differ inj-detail. 

With tbe-| rl^elp of rigid-iorcj and shell-model 
calculationgl3T3 some of the more pronounced structures 
can be connected with certain elongation patterns of the 
ions in the unit cell. For example, the peak in the /-In- 
curve at 4.64 THz in Tl2Ba2Ca2Cu30io seems to be due 
to a (Cu,Ba)-mode. 

The qualitative features of the subgap resonances have 
already been explained with the help of the phonon in- 
terpretation in Ref. ^ The position of the resonance is 
independent on temperature, magnetic field and the ge- 
ometry of the probe, while the intensity of the structure 
varies ~ with the critical current density jc{T,B). 
Also the behaviour of the position and intensity of the 
structures in external pressure are in agreement with the 
phonon explanation and formula (|l|). 

One of the main qualitative features of the general the- 
ory with dispersive phonon bands, which goes beyond the 
local oscillator model used in Refs. ^,||, is the possibility 
to describe resonances at van-Hove singularities, which 
appear e.g. at the upper band edge of the acoustical 
phonon band. 

This might be an explanation for a resonance seen in 
Ref. I at 3.2 mV (=1.54THz) in the /-y-characteristic 



of Tl2Ba2Ca2Cu30io, because the same frefluency is ex- 
pected by lattice dynamical calculationsE^ for the up- 
per edge of the acoustical band, and there are no optical 
phonon bands in this low frequency range. Fig. |^ shows 
a fit of a /-V-curve calculated with the two-atomic chain 
model (and some additional damping) to the experimen- 
tal results from Ref. 

Similarily in BSCCO the upper edge of the longitudinal 
acoustical phonon band, which hasJaeen detected in in- 
elastic neutron scattering at 5 me\o, might correspond 
to a less pronounced phonon resonance in the intrinsic 
Josephson effect at K^g = 4.9 mV (= 2.38 THz), which 
seems to be invisible in optical experiments. Neverthe- 
less at the present time this interpretation is not yet fully 
conclusive, as the instrumental resolution of 0.65 meV 
in the neutron scattering experiment is still rather large 
compared with the accuracy in the measurement of the 
I-V-cnvYQ. 
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FIG. 8. Fit of the experimental J-V-curve in 

Tl2Ba2Ca2Cu30io from Ref. ^ near the subgap resonances 
at the band edge of the acoustical branch (at 1.5 THz) and 
an optical branch with the help of the a two-atomic chain 
model. The inset shows experimental results over a wider 
frequency range. 

Also the effect of the two van-Hove singularities at the 
optical band edges on the /-T^-curve as discussed above 
might have been seen in the satellite structures at 5.17 
THz (10.7 mV) and 5.6 (11.6 mV) in the I-V-cnvYe of 
BSCCOEl. This interpretation is further supported by the 
theoretical prediction of the bandwidth ~ 0.3THz and 
the fact that the assignment of no other phonon mode is 
plausible. 
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TABLE I. Comparison of the frequencies /sg = ^Vdc (in 
THz) of the most pronounced subgap resonances in Ref. ^ 
and of infrared- and Raman active modes in Bi2Sr2CaCu2 08 
and Tl2Ba2Ca2Cu„02n+4. 



Subgap Resonances and Phonons 


in Bi2Sr2CaCu2 08 


/ag 2.97 3.89 5.17 5.60 Ref. 


M. 


Josephsoneffect 


/lo 2.85 5.07 Ref. 


11 


IR-reflecticity 


/lo 2.86 5.16 Ref. 


22 


IR-reflecticity 


/to 3.80 Refs. 




Raman effect 


Subgap Resonances and Phonons 


in 


TI2 Ba2 Ca2 Cu„ 02n+4 


/sg 3.63 4.64 Ref. 


M. 


Josephsoneffect^ 


/lo 4.50 Ref. |T 




IR-reflecticity'' 



"Tl2Ba2Ca2Cu30io 
''Tl2Ba2Ca2Cu208 



VII. CONCLUSIONS AND OUTLOOK 

In this paper we have developed a microscopic the- 
ory for the coupling between Josephson oscillations and 
phonons in intrinsic Josephson systems like the highly 
anisotropic cuprate superconductors. We determined the 
precise form of the longitudinal dielectric function Eq. 
( ^ ) describing this coupling and obtained analytical re- 
sults Eqs. ( [s^p^ ) for the I-V-cuive for one and several 
resistive junctions. The principle features and selection 
rules for phonon resonances in the I-V-cmves are illus- 
trated with help of a simple lattice dynamical model. 

We have shown that not only optical, but also acous- 
tical phonons at the edge of the Brillouin zone couple 
to Josephson oscillations. This may explain the struc- 
ture observed in Ref. ^ in the I-V-cwcve occurring at an 
unusual low voltage (frequency), which is not found in re- 
flectivity experiments, testing transverse optical phonons 
at g = 0, and in lattice-dynamical calculations for in- 
frared active phonons at g = 0. A weak,satellite-structure 
observed in the I-V-curve of BSCCOQ may be due to a 
double resonance from the two van-Hove singularites of 
an optical branch. 

The analytical results are compared with results from 
a numerical integration of the coupled equations of mo- 
tion for the Josephson oscillations and phonons. For this 
purpose a simplified lattice dynamical model has been 
used with one acoustical and one optical branch. It is 
found that in the limit of large values of the McCumber 
parameter the numerical results follow closely the ana- 
lytical solutions with the following exceptions: 1. Using 
a gradual change of the bias-current, regions of the /-In- 
curve with negative differential conductivity are skipped 
as is observed in current-biased experiments. 2. In the 
case of several resistive junctions, where several analyt- 
ical solutions are obtained, the numerical result follows 
only one of the analytical solutions. The stability of the 
different analytical solutions is currently investigated. It 
seems to be that the solution which gives a minimum 



for the interaction energy between polarisation and the 
electric field generated by the Josephson oscillations at a 
given frequency is most stable. The phonons thus serve 
to synchronise the Josephson oscillations in different re- 
sistive layers, which is important for the application of 
such systems as high-frequency devices. 

In this paper we have considered only the case of a cur- 
rent distribution in the c-direction which is homogeneous 
along the layers. We neglected all magnetic field effects 
assuming that all quantities are uniform along the lay- 
ers. For real systems this is an artificial approximation, 
because the current induced magnetic field can never be 
avoided completely. Nevertheless, we argue that this is 
a valid approximation for the intrinsic Josephson sys- 
tems forming a mesa structure of about fifty layers with a 
width of several fi and a thickness which is much smaller. 
Such a junction is still short with respect to the length 
Aj = c/ujp describing the variation of the phase-difference 
along the layers (in order to avoid completely the spon- 
taneous formation of vortices, which have a smaller di- 
ameter, the mesa width should not exceed 1-2 /i). There- 
fore the treatment of the superconducting layers as metal 
sheets with a uniform charge distribution, the creation 
of uniform polarisation fields and the neglect of g|| in the 
calculation of the longitudinal dielectric function is jus- 
tified for the systems investigated. This will be different 
in the case-pf longer junctions and strong external mag- 
netic fielder. In particular the flow of vortices and their 
intecaption with phonons has to be investigated in this 
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APPENDIX A: J-V"-CURVES FOR SEVERAL 
RESISTIVE JUNCTIONS 

Quite generally the RSJ-equation for the n-th junction 
can be written as 

j = jcsin7„(t) + CT^7n(0 (Al) 

t pOO 
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where epjj(n — n',t) is the dielectric function Eq. (31) in 
real space and time domain. 

Let us denote the index of a general junction by n G 
{1 . . . N} and the subset of indices of resistive junctions 
by /. In the following we assume that oscillations in all 
junctions are governed by one frequency uj and the phases 
can be written as 

l0n + 57nW else ^^^> 



with 



(57„(0 =<57„e-*"* + c.c. 



(A3) 



Expanding with respect to Sjn and keeping only the low- 
est harmonic one finds 



sin7„(t) 



sin(6lj + Lut) + ^{S-fje''>^ ) n = j el 
sm9n + cos6'„(57„(i) 



else. 



(A4) 



Splitting the RSJ-equation (Al) into a dc-part and an 
oscillating part one obtains for the dc-part: 

• ^ f crEdc + jc3?(<57je'^^ ) for n = j el x^^x 
\ jc sin 0n else, 

while the oscillating part can be written as (w^ = 

U^pSjnit) + -Sjnit) (A6) 



/oo 
e{n-n',t-t')S%,{t')dt' ^ Ut) 
-oo 



with 



/„ {t) = ( ^P^"fi (^) ~ + n^j el (^7) 

1 else. 



Writing 



/,(i) = /,e-'"* + , 



(A8) 



the amplitude of the driving term for a resistive junction 
is given by 



with the Green's function 

G-\q,,uj)^Ql-iuj^^~uj^el^^{q„Uj) (A12) 

of the homogeneous equation. 

From this an equation for {i e I) in the resistive 
junctions is obtained: 



2 



with 



Using 



(A13) 



(A14) 



e-ph(z,A:,c.) = -±G-\i,k,iu) + - —)S,,k (A15) 



we get 



eoLU 



i.k 



% (A16) 



which gives us the dc-component Eq. ( |A5| ) of the current 
density 



(A17) 



Jc 



fee/ 



Finally we want to note that the dielectric function 
eph(i, k,ijj) can also be used to write the RSJ equations 
in a form where only the phases 7i(i), of the resistive 
junctions i e I enter: 

j = jcsm-f,{t) + (^^yXt) (A18) 
+ ?4E r -^Ah3.t-t')^,{t')dt'. 



Introducing the spatial Fourier transform 
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1 ^ 



we obtain 

7fe) = G(<z„c.)^/,e-''^^^- 



(A9) 



(AlO) 



(All) 



^ R. Kleiner, F. Steinmeyer, G. Kunkel, P. Miiller, Phys. 
Rev. Lett. 68, 2394 (1992). 

^ R. Kleiner, P. Miiller, Phys. Rev. B 49 (1994), 1327. 

^ K. Schlenga, G. Hechtfischer, R. Kleiner, W. Walkenhorst, 
P. Miiller, Phys. Rev. Lett. 76, 4943 (1996). 

* A. Yurgens, D. Winkler, N. Zavaritsky, T. Claeson, Pro- 
ceedings of SPIE, Vol 2697, 433 (1996). 

'^P. Seidel, A. Pfuch, U. Hiibner, F. Schmidl, H. Schnei- 
dewind, T. Ecke, J. Scherbel, Physica C 293, 49 (1997). 



13 



^ Ch. Helm, Ch. Preis, F. Forsthofer, J. Keller, K. Schlenga, 

R. Kleiner, P. Miiller, Phys. Rev. Lett. 79, 737 (1997). 
^ Ch. Helm, Ch. Preis, F. Forsthofer, J. Keller, K. Schlenga, 

R. Kleiner, P. Miiller, Physica C 293, 60 (1997). 
* K. Schlenga, R. Kleiner, G. Hechtfischer, M. Mofile, S. 

Schmitt, P. Miiller, Ch. Helm, Ch. Preis, F. Forsthofer, J. 

Keller, H. L. Johnson, M. Veith, E. Steinbeifi, Phys. Rev. 

B 57, 14518 (1998). 
® T. Zetterer, M. Franz, J. Schiitzmann, W. Ose, H.H. Otto, 

K.F. Renk, Phys. Rev, B 41, 9499 (1990); V.M. Burlakov, 

S.V. Shulga, J. Keller, K.F. Renk, Physica C 203, 68 (1992) 

and references therein. 
^° R.G. Buckley, M.P. Staines, D.M. Pooke, T. Stoto, N.E. 

Flower, Physica C 248, 247 (1995) and references therein. 
" A. A. Tsvetkov, D. Dulic, D. van der Marel, A. Damas- 

celli, G. A. Kaljushnaia, J. I. Gorina, N. N. Senturina, B. 

Willemsen, N. N. Kolesnikov, Z. F. Ren, J. H. Wang, A. 

A. Menovsky, T. T. M. Palstra, Phys. Rev. B 60, 13196 

(1999) and references therein. 

P. Berberich, R. Buemann, H. Kinder, Phys. Rev. Lett. 49, 
1500 (1982). 

" Ya.G. Ponomarev, E.B. Tsokur, M.V. Sudakova, S.N. Tch- 
esnokov, M.A. Lorenz, M.A. Hein, G. Miiller, H. Piel, B.A. 
Animov, Solid State Comm. Ill, 513 (1999). 

" T. Koyama, M. Tachiki, Phys. Rev. B 54, 16183 (1996). 

^'^ Ch. Preis, Ch. Helm, J. Keller, A. Sergeev, R. Kleiner, 
SPIE Conference Proceedings "Superconducting Superlat- 
tices II: Native and Artificial", Vol 3480, 236 (1998). 
E.G. Maksimov, P.I. Arseev, N.S. Maslova, Solid State 
Comm. Ill, 391 (1999) . 

R. Zeyher, G. Zwicknagl, Solid State Commun. 66, 617 
(1988); R. Zeyher, G. Zwicknagl, Z. Phys. B- Condensed 
Matter 78, 175 (1990). 

A. D. Kulkarni, J. Prade, F. W. de Wette, W. Kress, U. 
Schroder, Phys. Rev. B 40, 2642 (1989). 
A. D. Kulkarni, F. W. de Wette, J. Prade, U. Schroder, W. 
Kress, Phys. Rev. B 43, 5451 (1991). 
2° J. Prade, A. D. Kulkarni, F. W. de Wette, U. Schroder, W. 
Kress, Phys. Rev. B 39, 2771 (1989). 
A. Yurgens, private communication. 

S. Tajima, G. D. Gu, S. Miyamoto, A. Odagawa, N. 

Koshizuka, Phys. Rev. B 48, 16164 (1993). 
23 C.-S. Jia, P.-Y. Lin, Y. Xiao, X.-W. Jiang, X.-Y. Gou, S. 

Huo, H. Li, Q.-B. Yang, Physica C 268, 41 (1996). 
^''H. A. Mook, B. C. Chakoumakos, M. MostoUer, A. T. 

Boothroyd, D. McK. Paul, Phys. Rev. Lett. 69, 2272 

(1992). 

Ch. Preis, Ch. Helm, K . Schmalzl, Ch. Wa lter, J. Keller, 



to appear in Physica C, cond-mat/0001150 
26 



' Ch. Preis, K. Schmalzl, Ch. Helm, J. Keller, in preparation. 



14 



